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Hemispheric  Asymmetries  Detected  From  Satellite 
Ionization  Gauge  Measurements 


1.  INTRODUCTION 

The  ionization  gauge  experiment  successfully  flown  onboard  the  Air  Force  S3-1 
satellite  resulted  in  the  accumulation  of  an  extensive  file  of  neutral  density  measure¬ 
ments  taken  over  the  period  extending  from  November  1974  to  May  1975.  This  re¬ 
port  is  a  first  attempt  at  performing  an  analysis  of  the  total  d  ita  file  gathered  as  a 
result  of  the  S3-1  flight.  The  satellite  orbit,  experiment  details,  reduction  proce¬ 
dures  including  computer  flow  charts,  and  representative  results  have  been  pre- 

1  2 

viously  presented  and  described  in  detail  in  earlier  reports.  *  The  statistical 
analysis  performed  here  is  similar  to  that  done  by  Forbes  et  al  in  which  the 
statistical  parameter  analyzed  is  the  ratio  R  between  the  measured  density  and  a 
model  density.  Frequency  distribution  plots  of  R  together  with  mean  values  and 
standard  deviations  are  given  for  different  levels  of  geomagnetic  activity.  The  three 
levels  of  geomagnetic  activity  used  are  characterized  as  quiet,  moderate,  and 
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disturbed  levels.  Comparisons  of  north  and  south  hemispheric  ratios  are  con¬ 
structed  over  a  range  of  different  altitudes  and  latitudes  for  the  levels  of  geomag¬ 
netic  activity  given  above.  The  ratio  analysis  provides  a  more  meaningful  com¬ 
parison  of  the  relative  hemispheric  behavior  of  atmospheric  density  than  would  a 
comparison  of  the  density  measurements  alone.  This  is  so  because  the  S3-1  satel¬ 
lite  measurements  were  obtained  over  relatively  short  spatial  intervals,  that  is, 
the  satellite's  data  acquisition  duty  cycle  routinely  involved  the  taking  of  one  perigee 
pass  every  alternate  revolution.  A  perigee  pass  consisted  of  a  20-min  tape  record¬ 
ing  centered  around  the  time  of  perigee  occurrance.  The  net  result  was  a  very  large 
accumulation  of  measurements  under  diverse  conditions  of  local  time,  latitude, 
longitude,  altitude,  season,  and  varied  geophysical  conditions  of  solar  flux  and  geo¬ 
magnetic  activity.  Figure  1  provides  a  time  history  of  the  satellite's  perigee  loca¬ 
tions  in  terms  of  altitude,  local  time,  and  latitude.  In  this  report  the  Jacchia  1971 
4 

model  is  used  to  form  the  ratio  parameters,  where  the  model  values  are  adjusted 
to  the  appropriate  conditions  pertaining  to  the  time  of  each  measurement.  An  ad¬ 
vantage  in  using  the  J71  model  for  forming  the  ratio  parameter  is  that  the  latitudinal 
variation  of  the  diurnal  bulge  in  the  Jacchia  model  is  symmetrical  in  both  hemis¬ 
pheres  wherein  the  maximum  density  migrates  north  and  south  of  the  equator  in 
phase  with  the  sub-solar  point.  A  disadvantage  to  the  use  of  any  model  in  deter¬ 
mining  the  ratio  parameter  is  that  whatever  model  is  employed,  inherent  in  the 
ratios  are  the  inadequacies  of  the  modelling  technique  as  well  as  unmodelled  real 
atmospheric  differences. 


2.  STATISTICAL  PROCEDURE 


The  three  parameters  derived  from  the  data  base  are  the  mean  ratio  R  (mean 
of  the  ratios,  data/J71),  the  standard  deviation  and  N  the  total  number  of  mea¬ 
surements.  The  expressions  used  to  determine  IT  and  o_  are  as  follows: 
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Observatory,  Cambridge,  AlA. 
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Figure  1.  Satellite  Position  in  Terms  of  Perigee  Location,  With  Local 
Time,  Altitude  and  Latitude  of  Perigee  Given  as  a  Function  of  Date  and 
Orbit  Number 


The  data  base  is  divided  into  different  combinations  of  parameters  to  form  data 

sets.  Within  each  set,  the  total  number  of  points  (ratios)  are  sorted  in  terms  of 

their  occurrence  and  partitioned  into  fixed  intervals  to  form  frequency  histograms. 

For  each  histogram  the  three  parameters  It,  (7^  and  N  are  computed  and  a  normal 

distribution  curve  constructed  which  is  then  superimposed  over  the  histogram.  The 

5 

expression  used  in  deriving  the  normal  distribution  curve  is: 

5.  Freeman,  H.  (1963)  Introduction  to  Statistical  Inference,  p  140,  Addison- 
Wesley  Publishing  Co. 
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where  Y,  the  ordinate,  is  the  number  of  occurrences  per  interval  i.  The  normal 
curve  superimposition  although  interesting  is  not  too  significant.  Its  worth  is 
greatest  when  used  for  model  evaluations.  Some  insight,  however,  is  gained  from 
histogram  examination  as  one  can  determine  how  results  are  distributed  about  their 
means  and  detect  asymmetric  distributions.  The  real  value  of  the  statistical  treat¬ 
ment  is  derived  from  the  employment  of  the  mean  ratios  and  the  variation  of  these 
ratios  under  diverse  combinations  of  solar,  geophysical  and  global  conditions 
obtainable  from  an  extensive  data  base.  Used  in  this  manner  the  statistical  analy¬ 
sis  provides  an  excellent  technique  for  investigating  average  atmospheric  behavior 
and  examining  systematic  variations. 


3.  DATA 

The  data  base  partitioning  used  in  this  report  is  shown  in  Figure  2  and  the 
respective  sample  size  for  each  bin  category  is  shown  in  Figure  3.  As  seen  in 
Figure  2,  the  base  is  first  separated  into  the  two  central  divisions  of  northern  and 
southern  hemispheric  data  followed  by  a  further  division  into  three  sub-levels. 

The  first  sub-level  separates  hemispheric  data  according  to  the  degree  of  geomag¬ 
netic  activity;  quiet,  moderate,  or  disturbed.  Or  more  precisely,  the  geomagnetic 
division  is  determined  by  the  value  of  the  Kp  planetary  index  (6  hr  lag).  The  range 
limits  for  the  activity  levels  used  are: 


Quiet  activity  (low); 

Moderate  activity  (moderate); 
Disturbed  activity  (high); 


0  £  K  -  3 

P 

3  <  K  ?  5 

P 

5  <  K  -•  9  . 

P 


The  second  and  third  sub-level  breakdowns  are  by  altitude  and  latitude  respectively. 
These  breakdowns  are  determined  in  a  straightforward  manner  reflecting  a  com¬ 
promise  between  adequate  statistical  sample  sizes  and  a  desire  to  avoid  an  excessive 
number  of  computer  runs  needed  to  compute  numerical  values  and  to  run  off  plotting 
routines.  The  partitioning  as  given  in  Figure  2  involved  the  generation  of  135  fre¬ 
quency  histograms,  a  portion  of  which  are  shown  in  the  appendix.  In  addition  to  the 
above,  various  bins  were  combined  to  provide  larger  bins  resulting  in  more  histo¬ 
gram  computations.  An  example  of  this  is  shown  in  Table  1  where  combinations  of 
several  bins  were  constructed  in  order  to  develop  a  more  meaningful  comparison 
between  gauge  results  shown  here  and  those  obtained  from  another  technique. 
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Figure  2.  Partition  Diagram  Outlining  Data  Base  Division 


Histograms  with  normal  curves  superimposed  are  included  in  the  appendix  together 
with  legend  data  consisting  of  values  of  71,  ar>,  X  and  the  normalized  V  value  as 
well  as  the  parameters  that  define  the  bin  confines.  Only  a  portion  of  the  total 
number  of  histograms  generated  as  a  result  of  the  bin  partitioning  has  been  included 
in  the  appendix.  However,  that  portion  is  representative  of  the  full  histogram  set. 
Figures,  which  are  discussed  in  detail  in  the  following  section,  representing  mean 
behavior  of  neutral  density  over  the  globe,  were  developed  that  show  results  as  a 
function  of  altitude  under  different  conditions  of  geomagnetic  activity  and  over 
several  altitude  regions. 
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Figure  3.  Sample  Sizes  Resulting  From  Data  Base  Partition  Given  in  Figure  2. 


4.  DISCUSSION 

Figure  4  shows  differences  between  mean  ratios  at  corresponding  latitudes  in 
the  north  and  south  hemispheres  plotted  as  a  function  of  latitude.  The  ratio  differ¬ 
ences  are  given  for  three  levels  of  geomagnetic  conditions;  quiet,  moderate,  and 
disturbed,  and  for  three  altitude  regions.  Among  the  prominent  features  seen  in 
this  set  of  curves,  perhaps  most  outstanding  is  the  total  absence  of  positive  differ¬ 
ences  below  200  km  over  all  latitudes.  Only  at  higher  altitudes  and  in  regions  about 
the  equator  do  we  observe  positive  ratio  differences.  The  shapes  of  the  difference 
curves  above  200  km,  also  appear  to  develop  into  similar  basic  patterns  that  tend 
to  be  accentuated  as  the  level  of  geomagnetic  increases.  This  structural  pattern 
is  manifested  in  greater  detail  in  the  later  figures. 


A{200-250km) 
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Figure  4.  Differences  in  Average  Density  Ratio  vs  Latitude  for  Three 
Geomagnetic  Groups.  Top  left  to  right:  low,  moderate  and  high  activity. 
Top  for  altitudes  less  than  200  km,  middle  for  altitudes  200  to  250  km 
and  bottom  for  altitudes  250  to  300  km 


The  overall  negative  differences  seen  in  Figure  4  indicate  that  atmospheric 
density  over  the  period  of  this  study  was  greater  over  the  southern  hemisphere, 
which  is  in  agreement  with  what  was  found  by  Barlier  et  al.  '  ^  His  examination  of 
over  12,000  values  obtained  from  drag  data  of  satellites  indicated  that  densities  in 
the  southern  hemisphere  were  systematically  larger  than  or  equal  to  the  corres¬ 
ponding  values  in  the  northern  hemisphere.  However,  in  terms  of  atomic  oxygen 


6.  Barlier,  F. ,  Perret,  D. ,  and  Jaeck,  C.  (1 973)  Structure  of  the  neutral  atmos¬ 
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North-south  asymmetries  in  the  thermosphere  during  the  last  maximum  of 
the  solar  cycle,  J.  Geophys.  Res.  ^5273-5285. 
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concentrations,  it  was  shown  that  at  200  km  the  concentrations  tended  to  accumu¬ 
late  over  the  cold  hemisphere  (winter)  at  the  lower  latitudes  where  seasonal  effects 
did  not  dominate,  which  also  is  not  in  disagreement  with  what  we  observe  in 
Figure  4. 

It  is  interesting  to  note  that  the  Jacchia  71  model  treats  hemispheric  variations 
in  a  uniform  manner,  that  is,  the  diurnal  maximum  migrates  symmetrically  about 
the  equator  in  phase  with  the  declination  angle.  In  forming  our  ratios  the  normali¬ 
zation  of  the  density  measurement  to  the  J71  model  value  should,  in  practice,  com¬ 
pensate  somewhat  for  the  seasonal  effects  due  to  the  migration  of  the  diurnal  bulge 
across  the  equator.  As  the  southern  hemisphere  measurements  in  our  data  base 
were  taken  almost  entirely  during  the  summer  season,  the  appropriate  southern 
hemispheric  J71  model  values  will,  all  other  factors  being  equal,  be  greater  than 
the  J71  northern  hemispheric  values  (spring/fall);  consequently,  from  the  model 
viewpoint  the  north  hemispheric  ratios  should  be  larger  than  south  hemispheric 
ratios. 

Figure  5  illustrates  global  latitudinal  structure  more  clearly;  shown  here  are 

ratios  of  measured  densities  to  model  densities  given  as  a  function  of  north  and 

south  latitudes  for  three  altitude  regions  and  different  levels  of  geomagnetic  activity. 

In  the  low  altitude  (<  200  km)  curve  set,  northern  ratios  are  fairly  constant  while 

southern  ratios  increase  in  the  direction  of  the  pole.  Again,  southern  ratios  are 

seen  to  be  considerably  larger  than  northern  ratios.  There  appears  to  be  little 

noticeable  effect  from  increasing  geomagnetic  activity  at  the  lower  altitudes.  For 

the  higher  altitude  curve  sets,  two  prominent  features  stand  out;  the  enhancement 

in  the  southern  auroral  region  during  periods  when  K  >  5  and  the  north-south 

P 

bulge/trough  structure  located  about  the  equator.  It  is  also  interesting  to  note 
that  for  the  most  part  the  ratio  magnitudes  appear  to  be  inversely  proportional  to 
the  level  of  geomagnetic  activity  indicating  the  J71  model  overcompensates  for 
geomagnetic  effects.  The  auroral  region  enhancement  is  present  only  in  the 
southern  hemisphere  auroral  region  where,  unlike  the  data  below  200  km,  the  den¬ 
sity  is  seen  to  decrease  poleward  of  the  auroral  region.  This  same  effect  was  ob- 

g 

served  by  Forbes  et  al  from  accelerometer  data  taken  onboard  the  NASA  AE-C 
satellite,  however,  in  that  report  north  and  south  latitudes  were  not  separated. 

This  particular  feature  found  here  and  in  the  AE-C  results  was  deduced  in  both 
cases  from  rather  small  statistical  data  samples,  therefore  the  observations  re¬ 
garding  the  auroral  region  feature;  enhancement  and  subsequent  poleward  drop-off, 

0 

should  be  viewed  as  inconclusive.  Rice  and  Sharp  have  reported  observations 

8.  Forbes,  J.  M, ,  Marcos,  F.A.,  and  Champion,  K.S.W.  (1977)  Lower  thermos¬ 

phere  response  to  geomagnetic  activity.  Space  Res.  XVIII:173-176. 

9.  Rice,  C.J. ,  and  Sharp,  L.  R.  (1977)  Neutral  Atmospheric  Waves  in  the 

Thermosphere  and  Tropospheric  Weather  System,  SaMSO-TR-^7-98. 


from  the  ion  gauge  measurements  performed  onboard  the  AE-C  spacecraft  wherein 
they  observed  an  asymmetric  auroral  region  feature  consisting  of  a  southern  lati¬ 
tude  region  (40°  to  68°S)  in  which  there  was  a  high  occurrence  of  wavelike  structural 
features  not  seen  in  the  corresponding  northern  latitude  region. 


o  LOW  Kp 
A  MODKp 
O  HIGH  Kp 


O  LOW  Kp 
A  MOD  Kp 
□  HIGH  Kp 


O  LOW  K„ 


A  M0DK- 
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GEOGRAPHIC  LATITUDE 

Figure  5.  Average  Density  Ratio  vs  Geographic  Latitude 
for  Three  Altitude  Regions  <  200  km.  200  to  250  km,  and 
250  to  3OT  km.  Open  circles  are  low  Kp  values,  triangles 
are  for  moderate  K  values  and  squared  are  for  high  K 
values  P  P 
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The  other  prominent  structure  observation  apparent  in  Figure  5  is  the 
equatorial  asymmetry,  a  north-south,  bulge /trough  formation  located  about  the 
equator.  This  bulge /trough  equatorial  feature  is  highly  visible  at  the  higher  alti¬ 
tudes  but  seems  to  disappear  below  200  km.  Certain  aspects  or  characteristics  of 
this  structure  are  more  evident  in  Figure  6  where  ratio  curves  are  grouped  to¬ 
gether  in  sets  according  to  the  level  of  geomagnetic  activity.  Foremost,  is  the 
phase  relationship  between  ratios  below  200  km  to  those  above  200  km.  The  equa¬ 
torial  structure  below  200  km  is  clearly  out  of  phase  with  the  structure  above 
200  km.  Below  200  km  at  northern  equatorial  latitudes  there  is  a  trough  in  place  of 
the  bulge  seen  at  altitudes  above  200  km.  In  a  like  manner,  the  deep  southern 
latitude  trough  at  high  altitudes  is  reversed  at  low  altitudes  and  although  one  does 
not  observe  a  well  defined  bulge  at  the  southern  latitudes  below  200  km,  one  can 
observe  a  shape  that  suggests  a  bulge  merged  with  other  effects.  A  second  charac¬ 
teristic  of  the  equatorial  structure  apparent  in  both  Figures  5  and  6  is  that  the 
minimum  of  the  southern  latitude  trough  appears  to  drift  towards  the  equator  during 

periods  of  increase  geomagnetic  activity  (K  >  5). 

P  in 

Equatorial  located  structure  of  a  similar  nature  has  been  reported  elsewhere  * 

and,  although  no  explanations  were  given  by  the  authors,  in  each  case  the  structure 
was  characterized  as  persistent  as  opposed  to  being  of  a  transient  nature.  The 
persistent  nature  of  the  structure  described  here  requires  a  non-transient  type 
mechanism  be  invoked  if  a  satisfactory  explanation  is  to  be  found.  Recall  that  the 
ratio  values  are  mean  values  representative  of  data  obtained  over  a  six  month 
period.  One  possible  explanation  would  involve  solar  tidal  variations  of  density. 

The  north-south  pattern  observed  in  our  data  set  is  suggestive  of  tidal  mode  propa¬ 
gation  into  the  thermosphere  with  antisymmetric,  semidiurnal  modes  dominant. 

1 2 

Semidiurnr.l  tidal  effects  have  been  theoretically  examined  by  Hong  and  Lindzen 
and  large  equatorial  latitudinal  variations,  particularly  at  solar  cycle  minimum, 
have  been  predicted.  Our  data  base  lacks  adequate  local  time  coverage  to  properly 
evaluate  tidal  oscillations  as  an  explanation.  However,  al'hough  not  shown  in  this 
report,  the  data  given  in  Figures  5  and  6  have  been  further  constrained  to  periods 
of  constant  local  time  (  ±  1  hr).  The  two  local  time  periods  were  at  0700  and 
2200  hours.  Under  this  added  constraint  of  constant  local  time  the  observed  struc¬ 
tural  features  remained  basically  unchanged  indicating  that  the  variations  in  the  data 
are  not  due  to  differences  in  local  time. 

10.  Philbrick,  C.  R. ,  and  Mclsaac,  J.  P.  (1972)  Measurements  of  atmospheric 

composition  near  400  km.  Space  Research  X\7I:743-749. 

11.  Ching,  B.K. ,  and  Carter,  V.L.  (1974)  Ion  gauge  measurements  of  latitudinal 

density  variations  at  night,  Geophys.  Res.  Lett.  ^93-96. 

12.  Hong,  S.  S. ,  and  Lindzen,  R.  S.  (1976)  Solar  semidiurnal  tide  in  the  thermos¬ 

phere,  J.  Atmos.  Sci.  33:135-153. 
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300  km 


Figure  6.  Average  Density  Ratio  vs  Latitude  for  Three  Levels 
of  Geomagnetic  Activity  as  Shown.  Altitude  symbols  are 
circles  <  200  km,  triangles  200  to  250  km  and  squares  250  to 
300  km 


A  point  that  requires  investigation  is  the  absence  of  the  structural  features 

13 

discussed  above  from  the  accelerometer  results  obtained  onboard  the  same  S3-1 
spacecraft.  In  an  attempt  to  ascertain  the  reason  for  this  discrepancy  the  gauge 
ratio  parameters  given  in  this  report  were  recomputed  to  coincide  with  the 


13.  Marcos,  F.A. ,  Mclnerney,  R.  E. ,  and  Fioretti,  R.  W.  (1978)  Variability  of 
the  Lower  Thermosphere  Determined  From  Satellite  Accelerometer  Data, 
AFGL-TR-78-0134,  AD  A‘058-982. 
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partitioning  used  by  Marcos  et  al  in  their  processing  of  accelerometer  results. 
Tnis  ne'.v  partitioning  of  the  gauge  data  involved  combining  the  two  altitude  bins 
of  <  200  km  and  200  to  250  km  to  form  a  single  bin  extending  from  perigee  to 
250  km  and  expanding  the  latitude  spans  by  15°  to  form  30°  bins.  Putting  results 
from  both  techniques  in  the  same  format  provides  a  means  for  directly  comparing 
results.  This  is  done  in  Table  1  where  results  from  both  data  sets  are  listed. 
Also,  the  data  shown  in  Table  1  are  plotted  as  a  function  of  latitude  (Figure  7). 
Referring  to  Figure  7,  one  readily  observes  that  the  figure  shows  two  pairs  of 
curves  for  two  levels  of  geomagnetic  activity  in  which  results  from  each  technique 
are  compared.  The  upper  and  lower  curve  sets  can  be  referred  to  as  representing 


high  and  low  geomagnetic,  respectively,  although  the  K  limits  for  both  sets  are 

13  P 

slightly  different.  Marcos  et  al  omitted  data  in  the  range  from  2.  6  =  =  3.  6 

in  order  to  emphasize  the  differences  between  the  high  and  low  data  groups. 

Also,  in  that  study,  there  was  no  separate  geomagnetic  activity  category  that  could 
be  characterized  as  moderate,  hence,  a  comparison  of  results  for  moderate  geo¬ 
magnetic  conditions  is  omitted.  Upon  examination  of  Figure  7  one  can  observe 
that  the  agreement  between  results  obtained  from  both  techniques  is  generally  good, 
particularly  in  the  equatorial  region  where  the  structure  reported  in  this  report 
has  been  observed.  The  region  of  greatest  discrepancy  is  in  the  high  southern 
latitudes  which  is  also  the  region  where  statistics  are  poorest.  Accelerometer 
results  also  indicate  higher  southern  hemispheric  densities.  One  fact  made  immedi¬ 
ately  evident  from  Figure  7  is  that  when  the  ionization  gauge  results  are  configured 
in  the  format  used  in  the  accelerometer  study,  the  results  reveal  little  of  the  struc¬ 
ture  reported  here;  hence,  indicating  no  general  disagreement  between  the  two  tech¬ 
niques.  Two  reasons  are  offered  for  the  lack  of  structure  seen  in  Figure  7,  namely, 
the  bin  limits  were  too  broad  (latitude  and  altitude)  and,  certainly  a  major  factor, 
the  altitude  dependent  phase  reversal  behavior  reported  here.  By  grouping  data 
from  below  200  km  with  data  above  200  km,  the  result  is  a  tendency  to  smear  or 
negate  the  structural  differences  between  the  two  altitude  regions. 
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Figure  7.  Ratios  vs  Latitude  for  Two  S3-1  Satellite  Experiments.  Triangles  are 
ionization  gauge  values,  circles  are  accelerometer  values 

5.  SUMMARY  AM)  CONCLUSION 

A  statistical  analysis  of  thermospheric  density  measurements,  in  excess  of 
150,000  values,  described  in  this  report  has  revealed  hemispheric  asymmetries. 

In  general,  at  least  for  the  period  of  this  study,  southern  hemispheric  densities 
are  found  to  be  greater  relative  to  the  Jacchia  71  model  than  northern  hemisphere 
densities.  Density  enhancements  are  observed  in  the  southern  auroral  region  under 
disturbed  geomagnetic  conditions  (K^  >  5)  that  are  not  seen  in  the  corresponding 
north  auroral  region.  A  non-transient,  equatorial-located,  structural  feature  is 
observed.  This  structure  which  is  asymmetric  about  the  equator,  is  seen  to  inten¬ 
sify  with  increasing  altitude  and  geomagnetic  activity.  A  tentative  explanation  is 
advanced  that  attributes  the  equatorial  behavior  to  atmospheric  tidal  effects. 

Further  analysis  of  these  results  is  planned  to  investigate  seasonal  influence  upon 
the  observed  equatorial  structure  and  determine  the  fine  spatial  resolution,  as  well. 


References 


1.  Mclsaac,  J.  P. ,  Champion,  K.  S.  W. ,  Mclnerney,  R.  E, ,  and  Delorey,  D. 

(1976)  Ionization  Gauge  Measurements  of  Atmospheric  Density  From  a  Low 
Altitude  Satellite,  AFGL-TR-76-0i  13,  AD  .4032373.  * 

2.  Mclsaac,  J.P. ,  Mclnerney,  R.  E. ,  and  Delorey,  D.  (1978)  Satellite  Ionization 

Gauge  Measurements  of  Atmospheric  Density,  AFG L-T R-78-026 1 , 

AD  A061613. 

3.  Forbes,  J.  M. ,  Marcos,  F.  A. ,  and  Gillette,  D.  F.  (1979)  An  Evaluation  of 

Thermospheric  Models.  AFGL-TR-79-0140,  AD  B36638. 

4.  Jacchia,  L.  G,  (1971)  Revised  Static  Models  of  the  Thermosphere  and  Exosphere 

With  Empirical  Temperature  Profiles,  Spec.  Rpt.  3 32,  Smithsonian  ~ 

Astrophys.  Observatory,  Cambridge,  MA. 

5.  Freeman,  H.  (1963)  Introduction  to  Statistical  Inference,  p  140,  Addison- 

Wesley  Publishing  Co. 

6.  Earlier,  F. ,  Perret,  D. ,  and  Jaeck,  C.  (1973)  Structure  of  the  neutral- 

atmosphere  between  150  and  500  km.  Space  Research  >(111:349-355. 

7.  Earlier,  F. ,  Bauer,  P. ,  Jaeck,  C. ,  Thuillier,  G. ,  and  Kockharts,  G.  (1974) 

North-south  asymmetries  in  the  thermosphere  during  the  last  maximum  of 
the  solar  cycle,  J.  Geophys.  Res.  |79:5273-5285. 

8.  Forbes,  J.  M. ,  Marcos,  F.  A. ,  and  Champion,  K.  S.  W.  (1977)  Lower 

thermosphere  response  to  geomagnetic  activity.  Space  Res.  XVIII:  173- 176. 

”  — 1 -  VVA/WW 

9.  Rice,  C.  J. ,  and  Sharp,  L.  R.  (1977)  Neutral  Atmospheric  Waves  in  the 

Thermosphex~e  and  Tropospheric  Weather  System,  sAmso-TR-77-98. 

10.  Philbrick,  C.  R. ,  and  Mclsaac,  J.P.  (1972)  Measurements  of  atmospheric 

composition  near  400  km.  Space  Research  XVl:743-749. 

11.  Ching,  E.  K. ,  and  Carter,  V.  L.  (1974)  Ion  gauge  measurements  of  latitudinal 

density  variations  at  night,  Geophys.  Res.  Lett.  1:93-96. 


References 


12.  Hong,  S.S.,  and  Lindzen,  R.S.  (1976)  Solar  semidiurnal  tide  in  the  thermos¬ 

phere,  _.L_Atmosi_Sci.  £3:135-153. 

13.  Marcos,  F.  A. ,  Mclnerney,  R.  E. ,  and  Fioretti,  R.  W.  (1978)  Variability  of 

the  Lower  Thermosphere  Determined  From  Satellite  Accelerometer  Data, 
AFGL-TR-78-0134,  AD  A058-982. 


Appendix  A 


Histograms  of  Frequency  Distributions  of  Ratios 
of  Measurements/Jacchia  71  Model 
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